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Using as a thermometer the temperature dependent magneto-transport of a two-dimensional elec- 
tron gas, we find that effective temperature scales with current as T a ~ I a , where a — 0.4 ± 2% in 
the Shubnikov-de Haas regime, and 0.53 ±2% in both the integer and fractional quantum Hall effect. 
This implies the phonon energy emission rate changes from the expected P ~ T 5 to P ~ T 4 . We ex- 
plain this, as well as the dramatic enhancement in phonon emission efficiency using a hydrodynamic 
model. 

PACS Numbers: 73.40. Hm, 73.50.Jt, 73.61.Ey 
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The concept of the inelastic scattering time is funda- 
mental to the physics of disordered quantum systems. 
Elastic scattering produces random interference patterns 
of the electron waves but does not introduce dccoher- 
ence. In contrast, inelastic scattering determines the 
phase coherence time which in turn controls the strength 
of the quantum interference observed for example in the 
Aharanov-Bohm effect in mesoscopic rings and in An- 
derson localization in disordered metals. Inelastic scat- 
tering is also fundamental to questions of energy equili- 
bration and cooling. In order to see Coulomb blockade 
and related effects in mesoscopic devices it is necessary 
to go to extremely low temperatures and to isolate the 
devices from their environment by means of large series 
resistances (R > h/e 2 ). || It turns out to be difficult 
to keep the electrons transiting these thin film resistors 
from falling out of equilibrium and heating up Q due to 
the long inelastic scattering time at low temperatures. |j| 

One difficulty in studying the dynamics of inelastic 
scattering in ordinary metallic films is the problem of 
making and comparing samples with significantly differ- 
ent characteristics. We overcome this problem here by 
studying high mobility two-dimensional electron gases 
both in the low magnetic field Shubnikov-de Haas (SdH) 
regime and the high field quantum Hall (IQHE and 
FQHE) regimes within the same sample. 

Previous models of cooling by phonon emission have al- 
ways assumed that the electrons are in plane wave states 
and invoke momentum conservation. Simple phase space 
considerations || then dictate that the power radiated 
into phonons is P ~ T^ d+2 \ where d is the dimension 
of space seen by the phonons (d = 3 in this case). The 
factor T d arises from the phonon phase space and the 
statically screened phonon matrix clement. [This applies 
both to metal films and inversion layers in piezoelectric 
media.] One additional factor of T represents the mean 
energy per phonon and the final factor of T represents 
the number of electrons which are sufficiently thermally 



excited to emit a phonon. Equating the radiated power 
to the Joule heating gives 
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I a with a = 0.40 ± 2% in 



Experimentally we find T c 
the SdH regime. 

In the presence of strong disorder or high magnetic 
field, the reduced mean free path (and magnetic length) 
can put the system into a new, hydrodynamic regime. 
Experimentally we find a new exponent a — 0.53 ± 2% 
(corresponding to P ~ T 4 ) and a nearly two order of 
magnitude enhancement of the emission rate. Mittal Q| 
has discussed the failure of existing theories to properly 
describe this regime. We present a new hydrodynamic 
theory of phonon emission in which dynamic screening 
plays a crucial role and which is in quantitative agree- 
ment with these unexpected experimental results. 

Current (I) or electric field (E) dependent transport 
measurements have been used extensively to study the 
electron phonon scattering process in 2 DEG's in separate 
regimes such as at zero B ||, in the SdH regime 0], in 
between Hall plateaus in the IQHE regime @Jll), and 
at the v = | FQHE plateau [l2| |. The latter indicated 
a possible dramatic enhancement in the electron-phonon 
coupling. However, no experiment has been performed 
that systematically explores different regimes in the same 
sample to explicitly investigate the consequences of the 
QHE, even though it has been suggested that the QHE 
will affect the electron phonon scattering [jl3| . 

The sample used in this work is a GaAs/AlGaAs 
heterostructure with an electron density n — 0.65 x 
10 11 cm" 2 , and a mobility \i ~ 500, 000 cm 2 /Vs a t 0.1K. 
The high mobility was chosen to allow the SdH oscilla- 
tions to be seen at low B. The donor spacer thickness 
is 1800A, and the Hall bar pattern on the sample has 
a channel width of 300 /im. The large device is used to 
avoid the influence of edge effects [|14| , which are known 
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to be more pronounced in small devices, jl5| and to avoid 
problems with heat being carried out of the sample di- 
rectly by electron diffusion. The sample and a ruthe- 
nium oxide resistor, to measure the temperature (T b ) at 
the sample position, are mounted in a dilution refrigera- 
tor [10, 11(c)] for 45mK< T b <1K. The transport coeffi- 
cients are measured by the standard ac lock-in technique. 
It is known that low-frequency ac and dc current produce 
the same results for the current dependence [10, 11(b)]. 

The inset of Fig.l shows p xx as a function of B 
for 0.06 T < B < 0.16 T (which corresponds to spin- 
unresolved Landau level filling factors, 18 < v < 38) at 
3 different excitation currents (I) at fixed T b = 100 mK. 
We measure the T b dependence of these SdH oscillations 
for 100 mK < T b < 750 mK while applying / = 5nA. 
By fitting the T dependent p xx amplitudes to Ando's 
semi-classical formula |l(J, we obtain the electron effec- 
tive mass 0.068 ± 0.001 in units of the free electron mass, 
comparable to that obtained from cyclotron resonance 
measurements |Hj]. 

Using the T dependent p xx amplitude as a thermome- 
ter, we find the same relation between T e and I for all 
oscillations down to 0.092 T in the inset of Fig.l. Oscilla- 
tions at lower B exist only in a very small range of I. We 
choose to plot, in Fig.l, T e vs. I for the two oscillations 
pointed to by arrows in the inset of Fig.l, because their 
amplitudes are measurable in a wider range of /. The 
closed and open symbols represent the data taken from 
the peak at B = 0.139 T and the dip at B = 0.133 T 
respectively. Different symbols represent data taken at 
different T&'s which are labeled by the side of the data set. 
For clarity, the open symbol data set has been offset by 
scaling by a factor of 1/3. For each T b , there is a current 
Io, below which T e remains constant, and above which 
T e merges into a single power law: T c ~ I a over about 
one and a half decades in I with a close to 0.4. The solid 
line under each data set has a slope of 0.4 and is drawn 
for reference. For comparison, we also draw a dashed line 
of slope 0.5. The experimental value of a is obtained by 
first collecting all data points which have T e > 2T b for 
each fixed T b , and then performing a linear least-squares 
fit to the resulting data points. We obtain a = 0.4 ± 2%, 
where 2% sets the statistical 68% confidence limit. 

We also deduce the absolute energy relaxation rate 
(l/"7m) by equating the input Joule heating to the cooling 
rate P ~ CT e /5ri n , where C is the specific heat of the 
2 DEC at B = 0. We find l/r in ~ 7.3 x 10 8 T 3 sec" 1 ^ 3 
at B = 0.13 T. Price has calculated l/r in = 2.92 x 
10 8 T 3 sec _1 K -3 for screened piezoelectric coupling in the 
absence of disorder. jl8| This agreement indicates that 
statically screened piezoelectric electron-phonon coupling 
controls the physics in the SdH regime. 

In the same sample, we also study the effect of / on p xy 
in the FQHE regime The inset of Fig. 2 shows p xy as 
function of B at 3 different Z's with fixed T b = 100 mK, 
between well developed FQHE states at v = ^ and |. 
There is no signature of higher order fractional states for 
i < v < | in our T range. We characterize the effect of 
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FIG. 1. Electron effective temperature T e vs. applied cur- 
rent I at T b = 100 and 200 mK, for B = 0.139 T (closed sym- 
bols) and B — 0.133 T (open symbols). The solid line under 
each data set has slope 0.4, the dashed line 0.5. Inset: p xx 
vs. B at T b = 100 mK for three different Fs. 

I by assigning a T e at each /, where T e is ob- 
tained by using (dp xy / dB) max as a thermometer. Here, 
(dp xy I dB) max is the maximum slope of p xy vs. B for 
| < v < | measured with I = 3nA. In Fig. 2, we plot 
T e vs. I in this region using open symbols, scaled by a 
factor of 1/6. Different symbols represent data taken at 
different T b s or upon different cooldowns of the sample. 
We also study the IQHE regime for Landau level index 
N = 2 i (5 < v < 6) and 2 | (4 < v < 5), where | and | 
represent spin direction. We do not study the IQHE tran- 
sition for N = 1 I and 1 f because there appears to have 
FQHE structures in p xy . The result for TV = 2 J. is plot- 
ted in Fig. 2 by closed symbols. There are no data points 
for I > 2 fiA for the FQHE data because (dp xy /dB) max 
is saturated to the classical Hall resistance. However, it 
is still clear that in both the IQHE and FQHE regimes 
T e behaves the same and there is a single power law, 
T c ~ J a , over more than one decade in / with a ~ 0.5. 
The solid line of slope 0.5 and the dashed line of slope 
0.4 are drawn for reference. We obtain the best fit values 
a = 0.53 ± 2% in the FQHE regime and a = 0.54 ± 2% 
and 0.53 ± 4% for N = 2 j and 2 | respectively. This is 
significantly different from that (0.40 ± 2%) in the SdH 
regime. Similar results have also been obtained in a sec- 
ond sample. This change in a can be viewed as a change 
in the cooling rate from P ~ T 5 in the SdH to P - T 4 
in the QHE. 

We note that in our sample the peak values of p xx are 
212, 210, 335, and 4620 Q/n for B = 0.139 T, N = 2 |, 
N = 2 t, and 1/3 < v < 2/5 respectively at T = 100 mK. 
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FIG. 2. Electron effective temperature T e vs. applied cur- 
rent / in the IQHE for 5 < v < 6 (closed symbols), and in 
the FQHE for 1/3 < v < 2/5 (open symbols). Each curve is 
labeled by the bath temperature. The dashed reference line 
under each data set has a slope of 0.5, and the solid line 0.4. 
Inset: p xy vs. B at Tb — 100 mK for three different 7's in the 
FQHE for 1/3 < v < 2/5. 

The p xx and hence the Joule heating at the fractional 
transition is 22 times that at the integer transition shown 
in Fig. 2, and yet the two curves are nearly identical. In 
addition, the current threshold in Fig. 2 for seeing heating 
at Th = 100 mK is about twice as large as for the SdH case 
in Fig.l. This indicates that, at the integer transition, 
cooling is about 4 times as efficient as in the SdH regime, 
and at the fractional transition, fully 80 times as efficient. 
One possible explanation for this dramatic enhancement 
and the P ~ T 4 power law is outlined in the following 
hydrodynamic model. 

The standard Fermi's Golden Rule expression for the 
rate of phonon emission |lj,^| can be rexpressed in the 
following form to yield the power radiated by the 2DEG 
at temperature T e into the lattice at Tb = 0: 

P = ^^ A (Q)|M A (Q)| 2 n B (/3 c ^ A (Q))G( <Z , WA (Q)), (2) 

A,Q 



where G(q,u>) is given by 
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c<j a (Q) is the phonon frequency for polarization A, q is 
the projection of the wavevector Q onto the plane of the 
2DEG, M\(Q) is the electron-phonon matrix element, 
iib is the Bose-Einstein factor, [3 C — 1/fcsTe, k = 12.9 
is the background semiconductor dielectric constant and 



e(q, ui) is the 2DEG dielectric function. The expression 
for G in terms of the conductivity <t xx is the standard one 
used in the theory of surface acoustic wave absorption. 


We assume that at high B and low T the system is in 
the hydrodynamic limit where typical phonon wave vec- 
tors and frequencies are negligible compared to charac- 
teristic scales over which the conductivity varies, so that 
(j xx is a real constant. This is justified from direct mea- 
surement |22) that at the IQHE critical point the peak 
value of a xx (at q ~ 0) is approximately independent of 
frequency up to at least 14 GHz (hu> ~ 0.7K). 

At low temperatures, piezoelectric coupling dominates 
HH so that |M A (Q)| 2 - l/Q. Using the fact that 
a xx is constant, simple power counting in Eq.(|J) yields 
P ~ T 4 . The physical interpretation of this is the fol- 
lowing. In the hydrodynamic limit, momentum conserva- 
tion is lost because of frequent collisions of the electrons 
with the disorder potential. The rate energy is absorbed 
from a phonon mode depends only on the square of its 
electric field and is independent of wave vector and fre- 
quency. Hence the 2DEG acts like a black (more pre- 
cisely a "gray") body and emits phonons with the usual 
T 4 spectrum. 

In the clean limit, the power counting is different since 
the conductivity is non-local: Re(l/a xx ) « q/qTF^KeoVF, 
where qtf is the Thomas-Fermi wave vector, wp is the 
Fermi velocity. This changes the temperature exponent 
to P ~ T 5 and correctly reproduces the usual Thomas- 
Fermi static screening result |18| for the prefactor. In 
the hydrodynamic limit, the static screening approxima- 
tion fails (since charge fluctuations are no longer high fre- 
quency plasmons at u ~ q 1 / 2 but rather relax slowly with 
z = 1 dynamics: lj — i(a X x/'2neo)q)- It is not the static 
compressibility of the electron gas that counts, but rather 
its limited ability to move charge quickly. The smaller the 
conductivity is, the poorer the screening is and hence the 
larger the emitted power. Using a xx = e 2 kFi c /h gives an 
elastic mean free path of £ e ss 1.9/zm at B = and so 
the system enters the hydrodynamic regime q£ < 1 only 



at temperatures below ~ 12mK. If however a xx ~ < 
the hydrodynamic regime extends up to about IK. 
Evaluation of Eq. (g) yields 



P = 



ha* 



-()"> / 10 2 I Watts/m^ 



(4) 



This expression, which contains no adjustable param- 
eters, yields the correct T dependence and (using the 
measured values of p xx and p xy ) lies below the exper- 
imentally measured power by only 17% (15%) for the 
IQHE (FQHE) transition. Considering the simplicity of 
the model, this level of agreement is quite good. Notice 
that one of the advantages of this formulation is that, 
unlike previous formulations, it is not necessary to make 
any assumptions about the unknown density of states in 
the interacting 2 DEG. 

Equating the Joule heating and radiated power we have 
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T c = 24.90w ra ) 1 / 4 (J/A/m) 1 / 2 K, 



which is ~ 4% below the observed values. This result 
explains why the curves for the fractional and integer 
transitions are nearly identical. The factor ((JxxPxx) 1 ^ 
(which is unity at B = 0) is a weak function of B in 
the sense that it is nearly equal for the two plateau tran- 
sitions studied. |23| Thus the increased Joule heating is 
precisely compensated by increased electron-phonon cou- 
pling due to reduced screening. That is, the threshold for 
heating occurs, to a good approximation, at fixed current 
rather than fixed power. 

We note that the result a = 0.53 ± 2% is close to the 
prediction a = 0.5 by Polyakov and Shklovskii in their 
theoretical model of T-scaling p4| . They consider the 
insulating IQHE plateau regime where Coulomb variable 
range hopping dominates the transport. The electron- 
phonon interaction is not directly relevant in that study. 

The present hydrodynamic heating model for the 
metallic critical point assumes that phonon emission is 
the rate limiting step and that the electrons are effec- 
tively in equilibrium at some temperature. An alterna- 
tive picture can also be developed in which one recognizes 
that there is a diverging correlation time at the critical 
point of the quantum Hall transition. |2^,^6| Assuming 
that it is this internal time scale rather than phonon emis- 
sion that controls the dynamics, it can be shown [^6| that 
a = z/(l + z) = 0.5, also in agreement with the present 
experiments (since for Coulomb interactions, one expects 
the dynamical exponent z = 1). More theoretical and ex- 
perimental work will be needed to distinguish these two 



models. [26| One caveat for the quantum critical picture 
is that, because of the high mobility, the sample does 
not show critical point power-law T-scaling in the linear 
response transport in the T range studied. 

Finally we note that in silicon deformation potential 
coupling yields P ~ T 6 and hence an exponent a = 1/3. 
This may well explain the apparent failure of dynamical 
scaling in the experiments of Kravchenko et al. [^7| on 
non-linear response in Si inversion layers where a = 1/3 
is observed. 
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